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We review here the tests of fundamental physics based on the dynamics of solar system objects. 



1 Introduction 



The solar system has the most well known dynamics with a very ancient history of observations. 
GR (GR) was first tested at the astronomical scale since the dynamical impact of GR has then 
astronomical effects. The Eddington first observations of the light bending during the 1919 
solar eclipse and the explanation by GR of the perihelia advance of Mercury in 1915 are the first 
stones of regular checks obtained at the astronomical scale with the Moon, the planets and the 
asteroids of the solar system. Since the nineties, planetary dynamics was drastically improved 
thanks to the numerous space missions orbiting planets and very accurately tracked by radio 
and VLBI observations. Since the Apollo missions, the Moon is also intensively observed with 
even greater accuracy. The figure Q] pictures the situation. For the Moon, the accuracy is at the 
centimeter level on a period of 35 years. For the planets, the geocentric distances are estimated 
at the meter level for Mars on a period of about 10 years due to MGS, MO and MEX tracking 
data. For Venus the improvement is also important but over a more limited period (2006-2010) 
thanks to the VEX mission. The improvement of the accuracy of Jupiter orbit stops with the 
latest available flyby data obtained with the Cassini mission in 2000. This latest gives very good 
constraints on the Saturn orbit which is now known with an accuracy of 20 meters over 3 years. 
GR is then confronted with very accurate observed positions of the Moon and the planets. The 
solar system is then an ideal laboratory for testing gravity. 

In the same time, theoretical developments ask to be tested in the solar system or forecast GR 
violations at the solar system scales. One can cite for example the violation of the equivalence 
principal by the standard models for unification of quantum physics and GR, the supplementary 
advances of planet perihelia and nodes expected by the MOND theories, the variation with 
distance of the PPN parameters (3 and 7 induced by dark energy P) or string theory P, ^j), 
variation of the gravitational constant G induced by dark ene rgy P) or scalar field theories p^) 
as well as supplementary accelerations due to dark matter P^), MOND theories P) or modified 
gravitational potentials P,^^). 
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Figure 1: Best fitted residuals for Moon and planets. The planetary residuals have been estimated with INPOPlOa 



' when the Moon residuals labelled (Muller 2010 ) have been estimated by^^ and those labelled (Williams et 
al. 2009) have been presented by^^ as well as the APOLLO median uncertainties. 



2 General overview 

2.1 Tests based on direct spacecraft dynamics 

The accuracy of the spacecraft tracking data orbiting a planet, in interplanetary phase or during 
a flyby can reach up to few centimeters over several months. Such tracking is done using doppler 
shift observations or VLBI differential measurements which are very sensitive to the deflection 
of light. With such accuracy, the tracking data of space missions seem to be a good tool to test 
gravity. However, some factors as navigation unknowns (AMDs, solar panel calibrations), planet 
unknowns (potential, rotation...), effect of the solar plasma, or the correlation with planetary 
ephemerides limit such gravity tests. Dedicated missions have then to be planed in order to 
overcome such difficulties. For example, the PPN 7 determination obtained by^^ was done 
during very specific periods of the Cassini mission, especially dedicated to such tests. 

The dynamics of the solar system planets and moons is less affected by non gravitational or 
unexpected accelerations and by technical unknowns and is constrained with also high accuracy. 



2.2 LLR tests 

With LLR observations, positions and velocities of the Moon are known with an accuracy from 
10 to 1 centimeter over 35 years. With the APOLLO project p^J), new developments in ob- 
servational techniques improve this accuracy with an observational error of about 1 millimeter. 
With such accuracy, plans improvements of at least one order of magnitude in the test of 
the equivalence principle, the determination of the PPN parameter (3 and of the test of inverse 
square law. The table Q] gathers the main tests of gravity done using LLR observations as well 
as planetary ephemerides and spacecraft tracking. The LLR analysis is clearly one of the main 
source of information for gravity. It produces tests with the best available accuracy for the 
equivalence principal, the prefered-frame tests and the detection of possible supplementary ac- 
celerations induced by dark matter. Present limitations in the modeling of the lunar interior, 
the Earth rotation as well as the planetary ephemerides induce differences between the teams 
analyzing the LLR observations ^33 E3) of several centimeters where the supposed accuracy of 
the APOLLO observations is about 1 millimeter. These discrepancies are obvious on figure [TJ 



Table 1: Results of gravity tests realized in the solar system. Columns 1 and 2 give the label of the tests and 
the studied objects. The Column 3 gives the obtained results with the mean value and 1 a least square standard 
deviations of the estimated parameter, except for the /3, 7, rb sup , and Cl sup obtained with the planets, see the 
text for more details. The last column gives the alternative theories of gravity wh ich can be constrained by these 
results. The values of /3 given here were all obtained with a value of 7 given by^l The ISL a obtained with LLR 



is for A = 4 x 10 km when the ISL a based on Mars data analysis is for A = 10 km. 
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An intensive work of comparisons and improvement of the Moon dynamics is now in progress. 



2.3 Tests based on planetary ephemerides 

As already discussed in section 1, the tracking observations of spacecrafts orbiting or flying by 
planets became very accurate in the nineties. The Mars data are crucial for the adjustment 
of the planetary ephemerides due to their high accuracy (about few meters on the Earth-Mars 
distances), their number (about 30 % of the complete data set) and their long time span (from 
1998 to 2010). This makes the planetary ephemerides very depend on this type of data and 
also very sensitive to the Mars orbit modeling and to the perturbations of the asteroids. On 
table [1] are found the tests done with planetary ephemerides. The dynamics of the planets give 
very competitive estimations of (5 and 7 PPN parameters as well as variations with time of the 
gravitational constant. Supplementary acceptable advances in the nodes and perihelia of planets 
are also constrained with the observations used in the adjustment of the planetary ephemerides. 
Tests of possible Pioneer-like accelerations on outer planet orbits have also been tested with 
high accuracy. The limitations of the gravity tests obtained with planetary ephemerides are 
mainly linked with the over weig ht of the Mars datasets and with the perturbations of the 
asteroids. As demonstrated kA^, there is a strong correlation on the geocentric Mars distances 
between the sun oblatness, the PPN parameter f3 and the mass of the asteroid ring used to 
average the perturbations of small asteroids on planet orbits. The decorrelation between the 
PPN parameters, the sun oblateness and the asteroid perturbations is then better obtained 
when the global adjustement of all the planet orbits is done simultaneously. Furthermore, the 
decorrelation between the /3 and 7 parameters are only possible if the two following equations 
are solved together in the fitting procedure. 

where Am is the advance of perihelia of the perturbed planet induced by GR (first term) 
and the oblateness of the sun J2 (second term with Rq, the sun radius). In the first equation, 
a,e are the semi-major axis and the eccentricity of the perturded planet and GMq is the mass 
of the sun. In the second equation, At is the Shapiro delay, the supplementary delay induced 
by the deflection of the light path by the sun. 



3 Results and discussions 

3.1 Equivalence principal 

A detailed description of the method used to test the equivalence principal in using LLR obser- 
vations is given in^. The test is an estimation of the gravitational mass to inertial mass ratio 
of the Moon and of the Earth. In GR, this ratio is equal to 1. However in PPN formalism, it 
can be formulated as 

Mo/Mr = l + ri[(^L)} 

where U is the body's gravitational self-energy, M is the mass of the body, c the speed of light 
and r\ a PPN parameter equal to in GR. In the equation of the geocentric motion of the Moon, 
0^ have introduced the differences in accelerations induced by Mq/Mj ^ 1 for the Moon and 
of the Earth. By comparisons to the LLR observations, it becomes possible to estimate the 
acceptable Mq/Mi ratio by direct least square fit of the numerically integrated acceleration or 
of the analytical estimations based on^4 The results presented in table Q] are the one obtained 
jjy[l5| combined with the laboratory estimations of the weak equivalence principal obtained by 
[22] 

As the PPN parameter 77 is linked to the (3 and the 7 by rj = 4/3 — 7 — 3, for a given value 
of 7, values for the /3 parameters can be deduced. Same methods have been applied by^with 
about the same results. Values obtained by^^ an are given in table [TJ 




Figure 2: Variations of postfit residuals obtained for different values of PPN /3 (x-axis ) an d 7 (y-axis). [1] stands 



for a PPN /3 value obtained by using LLR observations with 7 = 0, [2] stands for^ by a global fit of EPM 
planetar y ep hemerides. Kll stands for^^ determinations based mainly on Mars da ta analysis. M08 for^^ and 
W09 for^^ give values deduced from LLR for a fixed value of 7, B03 stands for^l determination of 7 by solar 



conjunction during the Cassini mission. 



3. 2 PPN parameter (3 and 7 

Since^ and^l, estimations of the PPN parameters are done with INPOP on a regular basis 
as well as estimations of acceptable supplementary advances of perihelia and nodes. A specific 
method presented irJ^is used for these estimations. In order to overcome the correlation prob- 
lems, we built up several planetary ephemerides for different values of the PPN parameters 
7) with a simultaneous fit of initial conditions of planets, mass of the Sun and asteroid densi- 
ties. On figure [2] are plotted the variations of postfit residuals induced by the modification of 
the corresponding (3 and 7 parameters. The left hand side plot gives the variations of postfit 
residuals including Mercury flyby normal point when the right hand side plot gives the varia- 
tions of residuals without the Mercury observations. The different levels of colors indicate the 
percentage of variations of the postfit residuals compared to those obtained with INPOPlOa. 
By projecting the 5% area on the /3-axis (or the 7-axis), one can deduced the corresponding (3 
(or 7) interval given in tabled] in which the residuals are modified by less than 5%. In looking 
at the two figures, one can see that the use of the Mercury flyby data give smallest intervals 
of possible /9,7- This is consistent with the fact that the Mercury observations are far more 
sensitive to gravity modifications than other data (see table 1 iir^ ) . 



3.3 Frame- dragging and prefered frame tests 

Based on LLR observations,^^ estimate a supplementary advance in the node of the Moon 
induced in GR by the motion of the Moon in the gravitational field of the Earth. The effect is 
called the de Sitter effect and the results are presented in table [U Prefered-frame coefficients 
a\ and ai have also been estimated by the same authors. Results are presented in table [TJ 
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Figure 3: Postfit residuals in meters obtained by comparison between Cassini tracking data and several 
ephemerides with added Pioneer-like accelerations. The characteristics of the accelerations are given in m.s -2 



3.4 Supplementary accelerations 

Supplementary accelerations can be induced by dark matter and dark energy surrounding the 
solar system or inside the solar system. Alternative descriptions of gravity can also induced 
modifications in the gravitational potential and then supplementary acceleration in the motion 
of solar system bodies and spacecraft. Possible tests have been made in introducing either 
constant accelerations in one specific direction (Pioneer-like anomaly with outer planets, dark 
matter with the Earth-Moon orbit) either accelerations induced by f(r) gravity or exponential 
potentials (ISL) . H51 an( j[IZl have constrained the ISL potential for the geocentric Moon (A = 
4 x 10 8 km) and Mars (A = 10 10 km). Some other estimations should be investigate for A > 10 12 
km. Results are given in table[TJ For the Pioneer-like accelerations, the figure [3741 gives the postfit 
residuals obtained by comparisons between the observed Earth-Saturn distances deduced from 
the Cassini tracking data and planetary ephemerides integrated with supplementary constant 
accelerations and fitted over the INPOPlOa data sample. These accelerations are similar in 
direction to the Pioneer accelerations but their modules vary from 10 -11 to 5.10~ 13 m.s- 2 . As 
one can see on figure I3T41 only accelerations smaller than 5.10 -13 m.s -2 are acceptable compared 
to the present accuracy of the Cassini observations. This result was confirmed by^-l 



3.5 Supplementary advances of perihelia and nodes 

New theoretical models [SIZS) forecast supplementary advances in the orbits of the solar system 
objects. Supplementary advances in the perihelia and nodes of the planets (from Mercury to 
Saturn) have also been tested in INPOP and are presented on table [2j All these results are 
based on the method presented in^l With INPOPlOa^H no supplementary advances in node 
or perihelia is detected when with INPOP08^3I supplementary advance in the Saturn perihelia 
was possible. Such variations can be explained by the improvement of INPOPlOa outer planet 
orbits compared to INPOP08. 



Table 2: Values of parameters obtained in the fit of INPOP08^31 and INPOPlOa to observations. The 
supplementary advances of perihelia and nodes are estimated in INPOPlOa and INP OP08 as the inter val in 
which the differences of postfit residuals from INPOPlOa are below 5%. P09 stands for P^) and P10 for p21). 
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4 Conclusions 

With the present gravity tests done in the solar system, GR is confirmed at the 10 -4 accuracy for 
PPN parameter (3 and the equivalence principal and at 10 -5 for PPN 7. No supplementary ad- 
vances of perihelia and nodes are detected at the present accuracy of the planetary ephemerides. 
Variations of the gravitational constant are not expected above 10 -13 yr _1 and stringent limits 
for the ISL tests are given for the inner solar system. Messenger tracking data would bring 
important informations for PPN parameter determinations and ISL tests should be done in the 
outer solar system. 
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